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Abstract 

The discrete Chebyshev polynomials t n (x,N) are orthogonal with respect to a 
distribution function, which is a step function with jumps one unit at the points 
x = 0, 1, • • • , N — 1, N being a fixed positive integer. By using a double integral 
representation, we derive two asymptotic expansions for t n (aN, N + 1) in the double 
scaling limit, namely, N — > oo and n/N — > b, where b G (0,1) and a G (—00,00). 
One expansion involves the confluent hypergeometric function and holds uniformly 
for a G [0, i], and the other involves the Gamma function and holds uniformly for 
a G (—oo,0). Both intervals of validity of these two expansions can be extended 
slightly to include a neighborhood of the origin. Asymptotic expansions for a > i 
can be obtained via a symmetry relation of t n (aN,N + 1) with respect to a = ^. 
Asymptotic formulas for small and large zeros of t n (x, N + 1) are also given. 
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1 INTRODUCTION 

In Szego [U p. 33], the discrete Chebyshev polynomials are defined by 
(1.1) t n (x, N) = n\A n ( X )( X ~ N ), n = 0, 1, 2, • • •, N - 1, 



K nj \ n 

where A denotes the difference operator 



;i.2) 



Af(x) = f(x + l)-f(x), 
A n f(x) = A{A n - 1 f(x)} 



n 



f(x + n)- ( 1 )f(x + n-l) + --- + (-l) n f(x). 



These polynomials are orthogonal with respect to the distribution da(x) of Stieltjes type, 
where a(x) is a step function with jumps one unit at the points x = 0, 1, • • •, N — 1, N 
being a fixed positive integer. More precisely, we have 



(1.3) / t n (x,N)t m (x,N)da(x)= t n (x,N)t m (x,N) =0 

J 00 x=0,l,-,N-l 

if n 7^ m, and 

/oo 
{t n (x, N)} 2 da(x) = £ {t n (x, N)} 2 
-OO n 1 !\r i 



;i-4) 



x=0,l,-,N-l 



N(N 2 - 1 2 )(N 2 - 2 2 ) ■ ■ ■ (N 2 - n 2 ) 
2n + 1 



for n, m = 0, 1, 2, • • •, N—l. The last two formulas hold for all non-negative values of n and 
m; in fact, they are trivial for n > N or m > N, since t n (x, N) = for x = 0, 1, 2, N—l, 
if n > N. 

The discrete Chebyshev polynomials can also be defined as a special case of Hahn 
polynomials 

Q n (x; a, (3, N) = 3 F 2 (-n, — x, n + a + /3 + 1; —AT, a + 1; 1) 

(1.5) = y> (-n) k (-x) k (n + a + /3 + l) k 

h> (-N) k (a + l) k k\ 

and we have 

t n {x, N) = (-l) n (N - n) n Q n {x- 0, 0, N - 1) 

(1.6) =(-in N -n) n <r { ~ n)k{ ~ x)k{n + 1)k - 



k=0 



(-N+l) k k\k\ 



see [2l p. 174 and 176] and [3]. These polynomials are used in least squares polynomial 
approximation jH p. 350] and sequential smoothing of numerical data [5]. (In jl], they are 
referred to as Gram polynomials; see Chihara [3], p. 162].) 
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In this paper, we are concerned with the problem of finding asymptotic expansions for 
the polynomials t n (x,N). Already, a tremendous amount of research has been carried 
out on the asymptotics of various discrete orthogonal polynomials. For instance, we have 
[0] and [7] for Charlier polynomials, [5] and [H] for Meixner polynomials, [TO] and [UJ 
for Krawtchouk polynomials, and [12] and [13] for Tricomi-Carlitz polynomials. There 
is also an important research monograph by Baik et al. [H] on a general asymptotic 
method, based on the Riemann-Hilbert approach, which in principle can be applied to 
all discrete orthogonal polynomials. However, their result is so general that it does not 
provide explicit asymptotic formulas for specific polynomials. Strangely enough, for the 
discrete Chebyshev polynomials there does not seem to be any asymptotic result in the 
literature; at least, we have not been able to find any, and the purpose of this paper is to 
fill in this gap. 

For convenience, we shall consider the behavior of t n (x,N + 1) instead of t n (x,N). 
Since the Hahn polynomials satisfy the equation [T5], (1.15)] 

fi V\ n (at a jva Qn{x;a,P,N) 

(L7) Qn ( N - X ' a >P> N) = Q n (N;a,P,N) 

and Q n (N; 0, 0, N) = (— l) n (note that there are slight differences in the definition of Q n 
[T5| (1.1)]), the discrete Chebyshev polynomials enjoy the symmetry relation 



(1.8) t n (x, N+l) = (-l)X(iV -x,N + 1). 

In view of (|1.8p . we may restrict x to the interval — oo < x < hN. We are interested in 
the behavior of t n (x, N + 1) as n, N — > oo in such a way that the ratios 

(1.9) a = x/N and b = n/N 
satisfy the inequalities 

(1.10) - oo < a < ~ and < b < 1. 

As a special case of our general result, we will show that for fixed value of x G [0, oo), 
there exists a number r\ < such that 

(-l)n+i T r n + N + 2 )N x n- 2x - 2 T(x + 1) 

am ^• Ar+1 '= F(Fttk 



x < sin TTX 



;+°(i)] +0 (^)}. 

This result also holds when x is negative, except that the term 0{e r ' N ) is absent. 



The arrangement of the present paper is as follows. In £j2j we present two (2-dimensional) 
integral representations for t n (x, N + l), one for the case x < and the other for x > 0. In 
§3j we transform the integration variables (w,t) to two new variables (u,r), and deduce 
the double integrals to their canonical forms, from which uniform asymptotic expansions 
can be derived by a standard integration-by-parts procedure. In §H we prove that the 
transformation (w, t) — > (u, r) is one-to-one and analytic. In §5], we investigate the analyt- 
icity of the integrands of the double integrals in (13.201) and (13. 29ft . The desired asymptotic 
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expansions are then derived in $6], and the estimation of the error terms is done in In 
§H1 we state the final results and give a list of special cases; and we extend the regions 
of validity of the final results in In ^TU| we present some asymptotic formulas for the 
large and small zeros of the polynomial under consideration. 



2 INTEGRAL REPRESENTATION 

From (jl.6p . we have 

t n (x, N + 1) = (-l) n (N + 1 - n) n Q n (x; 0, 0, N) 
(2-1) = { _ 1)n{N + 1 _ n)n J2 (Z?W-^ + !)* 



(-N) k klk\ 



Using beta integral, we also have 



n + l) k _ f _ 1)k r(n + 1 + k)T(N + 1 - k) 



, 22) (-N)k v ' r(n + i)r(iv + i 

. k T(n + N + 2 



-1 



r(n + 1)F(N + 1) 



/ (1 - t) n+k t N ~ k dt 
Jo 



for < k < n < N. Substituting ( 12. 2ft in ( 12. ID and interchanging the summation and 
integration signs, we obtain 

x , T(n + N + 2) 

ux,N+D=(-ir- 



( . ):} , T(n + l)T(N + l-n) 

(l-t) n t N 2 F 1 (-n,-x;l;l-r 1 )dt 

In the above integral, we now use the new scales a = x/N and b = n/N introduced in 
(jl.9p . and consider separate cases (i) < a < | and (ii) a < 0; cf. (jl.lOp . 

Recall the identity [16, (15.8.1)] 

(2.4) 2 F 1 (-n, -x- 1; 1 - r 1 ) = r^i^-n, 1 + x; 1; 1 - t) 
and the contour integral representation [2J p. 287] 

(2.5) 2 F 1 (a,&;c;z) = r(c) ^ r ^ ) ~ C) ^ ^^-^^(l-w)^, 

where Re 6 > and the curve 71 starts at w = 0, runs along the lower edge of the positive 
real line towards w — 1, encircles the point if = 1 in the counterclockwise direction, and 
returns to the origin along the upper edge of the positive real line. For case (i) < a < |, 
we apply ( 12 .4p and (12. 5p to (12. 3p . and obtain the double integral representation 

(2 6) t (x N + 1) = (-1)" r(n + iV + 2) A _J_ e ^ )(M 

1 J " l ' + j 27T* r(n + l)r(iV - n + 1) 7 7 7l w -l e 
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where 

(2.7) f(t,w) = 61n(l - t) + (1 - b) hit + alnw - aln(w - 1) +Mn[l - (1 - t)w). 

The function f(t,w) has a cut (— oo, 1] in the u>-plane. Since the coefficients of Int, 
ln(l — t) and ln[l — (1 — t)w] in Nf(t, w) are all nonnegative integers, there is no need to 
have cuts in the t-plane for Nf(t, w). 

For case (ii), we use instead of (I2.5P the contour integral representation 

(2.8) 2 F 1 (g, b- c; z) = ~^ c (1 _~ } &) / (-^^(l - w) M {l - wz)- a dw, 

where Re c > Re b and the curve 72 starts at w — 1, moves along the upper edge of the 
positive real line towards w = 0, encircles the origin in the counterclockwise direction, and 
returns to w — 1 along the lower edge of the positive real line. (A proof of this formula 
can be given along the same lines as that of (12.51) .) Thus, we have 

(2.9) t n {x, N + l) = Izff _^±^) f f 1 e *» )dmdt , 

2%i T(n + l)T(N - n + 1) J J 12 w - 1 

where 

(2.10) f{t,w) = 61n(l-t) + (l-b)]nt + aln(-w) -aln(l-w) + 61n[l - (l-t)w\. 

Similar to f(t,w) in (12. 7p . the function f(t,w) has a cut [0, +oo) in the w-plane and no 
cut in the t-plane. 



3 REDUCTION TO CANONICAL INTEGRALS 

To obtain large-A behavior of higher dimensional integral 

/(A) = [ f(z)e xs ^dz, z = (z u z 2 , ■ ■ ; z n ) e C\ 

where dz = dzi ■ ■ ■ dz n and 7 is a smooth manifold of (real) dimension n, Fdeoryuk [17] 
showed that just like the classical method of steepest descent, the saddle points of the 
phase function S(z) (i.e., the zeros of the first-order derivatives of S(z)) play an important 
role; see also Kaminski [18J. So, let us first investigate the location of the saddle points 
of the phase functions f(t,w) and f(t,w) in (12.71) and (12.101) . 

Upon solving the equations 

df df 
^(t,w) = and ^(^) = ° 



and the equations 



df df 

— (t,w)=0 and —{t,w) = 0, 
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one finds that the phase functions f(t,w) and f(t,w) have the same two sets of saddle 
points 



(3.1 
and 
(3.2 



2 - 2a - b 2 + bVb 2 - 4a + 4a 2 
2(l-a)(l + 6) : 

2 - 2a - b 2 - bVb 2 - 4a + 4a 2 



b + y/b 2 -4a + 4a 2 
25 



b - Vb 2 - Aa + 4a 2 



2(1 - a) (1 + 6) ' " 26 

Note that w + and to_, as well as £+ and coalesce when a approaches a + and a_, where 



(3.3) a± 
The following facts are easily verified: 



1. For a < 0, w± are both real and W- < < 1 < w + . As a — > 0, W- — > and 
iu + -)■ 1 + . 

2. For < a < a_, both w + and are real and < w_ < \ < w + < 1. As a — > a_, 
both w+ and approach |. 

3. For a_ < a < a + , w-t are complex conjugates, the real parts of which are both equal 
to At a — |, w-t attain their highest and lowest point on the line Re w = |. 

4. The locations of w± in the cases a + < a < 1 and a > 1 are similar to those in the 
cases < a < a_ and a < 0, respectively. This fact also reflects in the symmetry 
property of the discrete Chebyshev polynomial t n (aN,N + 1) in a with respect to 
a = |; cf. dUBD and ffTTOj) . 

Figure [1] summarizes the movements of w± . 
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Figure 1: Movements of to±. 
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In view of (11.101) . we only need consider the case — oo < a < |. We again divide our 
discussion into two subcases: (i) < a < | and (ii) — oo < a < 0. 

Case (i) < a < \. 

Returning to the integral representation (12. 6p . we now show that this double integral 
can be reduced to a canonical form, and that the original integration surface [0, 1] x 71 
in the t x w plane can be deformed into a new integration surface. An outline of the 
procedure is as follows. First, for each fix w G 71, we find a steepest descent path of 
the variable t, which passes through a relevant saddle point t (w) depending on w. We 
denote this path by p t . Note that not only the saddle point to(w), but all points t on the 
path p t depend on w. Thus, the phase function f(t (w), w) is now a function of w alone. 
Second, deform the integration path 71 in (12.61) into a steepest descent path of w, passing 
through w± given in (13.11) and (13.21) . We denote this path by "f w . It will be shown that 
these saddle points are indeed saddle points of f(t (w),w). 



To find the relevant saddle point to(w), we solve the equation 



and obtain 

(3.4) /-;ur 



± 2w - 1 ± y/l + Ab 2 (w - l)w 



2(1 + b)w 



Here, we have introduced two more cuts in the w-plane, both on the vertical line Re w = ~, 
one from (b+iy/l — b 2 )/2b all the way to w = |+zoo, and the other from (6— i\/l — b 2 )/2b 
all the way to w = ^ — zoo. As remarked earlier, the two saddle points (t + ,w + ) and 
(t_, W-) in (13. ip and (13. 2p play an important role in determining the asymptotic behavior 
of t n (x,N + 1). So, we must be careful in choosing the right saddle point from the two 
given m (ESP - For a < |, it can be shown that tQ(w±) = t±. Thus, t^w) is our relevant 
saddle point. Throughout the following discussion, we shall simply write t (w) for (w), 
when no misunderstanding would arise. 



As we have explained previously, for fixed w in the whole w-plane and b G (0, 1), the 
function e 1 ^^ 1 '^ inside the integral in (12. 6p has no singularity in the t-plane. However, 
to make it possible for us to deform the original integration path [0, 1] into a steepest 
descent path p t , we have to define the function f(t, w) in such a way that it avoids the 
cut (—00, 1] in the w-plane, and three additional cuts in the t-plane, (—00, 0], [1, 00) and 
also the radial line starting from 1 — 1/w and satisfying 1 — (1 — t)w < 0. With this in 
mind, we now define the steepest descent path p t , passing through t (w), by 

, . Im /(*, w) = Im f(t (w), w), 

Re f(t,w) <Re f{t {w),w); 

see Figure [2] and Figure |3j It is now standard to define a transformation t — > t by 

(3.6) r 2 = f(t (w),w)-f(t,w). 
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Figure 2: Steepest descent path in 
the t-plane when w is real. 




Figure 3: Steepest descent path in the 
t-plane when w is complex; e.g., w = 
0.5 + 0.2i. 



Note that we have r = — oo when t = 0, and r = +oo when t — 1. Furthermore, this 
mapping is one-to-one with r running from — oo to oo, and takes t = to(w) to r = 0. 
From (12. 7p and ( 13 . 61) . we have 



(3.7) 



dt 
dr 

dt 
dr 



2r 



2rt(l - t)[l - (l-t)w\ 



-df(t, w)/dt (1 + b)w{t - t+(w)){t - t (w)) ' 

N 1/2 

2t+H(i-t+H)[i-(i-t+HH 



Coupling (12. 6p and f 13 . 6 jl gives 

+ F( " + iV + 2) 

Z7T?, I in 

(3.8) 



0. 



27ri r(n+l)r(JV-n + l) 

+oo f y 



X 



oo »/ 71 



dt 



e 7V/( t0 ( W ),») e -^- dwdT) 

w — 1 dr 



where 71 is the integration path of w in (12.61) . Note that the first variable of / in 
(13. 8 j) is now to(w), instead of t. So the phase function depends only on w. Setting 
d f (t (w) , w) / dw = 0, we obtain saddle points 



(3.9) 



w± 



b±Vb 2 -4:a + Aa? 

2b : 



cf . (13. ip and (13.21) . We now deform 71 into steepest descent paths passing through w + and 
Recall that the movements of w± are summerized in Figured! as a runs along the 
real line. In the case a_ < a < |, the saddle points w± are complex conjugate numbers 
with real part equal to |. The steepest descent path passing through w + is given by 



(3.10) 



Im f(t (w),w) = Im f(t {w+),w + ), 
Re f(t (w),w) < Re f(t (w+),w + ). 



The steepest descent path through W- can be obtained by using the symmetry property 
with respect to the real axis; see Figure 0J 
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When < a < a_, both w± are real and lie in the interval [0, 1]; thus, w± appear on 
both upper and lower edges of the cut along the real line. We denote these points by 
u4 + ^ and w± ; see Figure |5j Since Re f(t (w_), w_) > Re f(t (w+ ), w+ ), w_ is a 
local maximum and the steepest descent path through w_ is along the real axis. From 
Figure [U it is clear that this path is on the real interval (0,w+^). The steepest descent 
path throught is the curve perpendicular to the real axis; cf. Figure EJ 

Motivated by the work of Jin and Wong j8], the movements of the saddle points w± 
and the steepest descent paths passing through them suggest that we should compare the 
integral of w in (13. 8p with the confluent hypergeometric function [161 P-326, (13.4.9)] 

(3.11) M(d, c, z) = EMKl+^Lg) f U d-H u _ i)^- V u dw, 

where Re d > and the integration curve 71 is the same one as that given in (12.51) or (13.81) . 
Since M(d, c, z) is a meromorphic function of c with poles at 0, —1, —2,- ■ -, the function 
M(d, c, z) := M(d, c, z)/T(c) is entire in both d and c; see [HI p.255]. With d = aN + 1 
(0 < a < |), c= 1, z = r}N, §FHi) gives 



(3.12) M(aN + 1,1, r]N) = ^- [ e 



71 U ~ 1 ' 



where 

(3.13) V , ( M ) = a In m — a ln(w — 1) + r/u 
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4. If, 



Figure 6: Steepest descent path of u when a > and —4a < rj < 0. 



)..( 




Figure 7: Steepest descent path of w when a > and r] < —4a. 



and the saddle points of ip are given by 



(3.14) 



u± 



rj ± -^/rf + 4ar] 
2rj ' 



The shape of the steepest descent paths of ip depends on the sign of a and the value of 
7]. Here, we are concerned with only the two cases: (i) a > and —4a < rj < 0, (ii) a > 
and rj < —4a; see Figures [6] and [0 

We also observe that when a > 0, the movements of u± are like those of w T , with the 
corresponding subcases 



a_ < a < - <H- -4a < ri < 0, 
0<a<a_ -H- 77 < —4a. 

Thus, we define the mapping w — > u by 



(3.15) 

with 
(3.16) 



f(t {w),w) = ip(u) +7 

= alnw — aln(w — 1) + 77W + 7 



u{wJ) = u + , 



where 7] and 7 are real numbers. Note that since w± are on both edges of the real line 
(similarly, u± are also on both edges of the real line) when < a < a_, (13. 16ft should be 
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understood in the sense: 

(3.17a) = u M , 

(3.17b) u{w { + } ) = u { S ] , 

It can be easily verified that (13 . 1 7aj) and (13.17bl) are equivalent in the case < a < a_. 
The existence and uniqueness of the numbers r] and 7 as well as the one-to-one and ana- 
lytic property of the mapping w — >■ u defined in (13.151) will be established in the following 
section. 



U[W 



(-h 



u 



(+) 

+ ' 

(-) 



Note that we have chosen the coefficients of lnw and ln(u — 1) in ip(u) to be a. The 
reason is that with this choice, the local expansions of f(t (w),w) at w = and w = 1 
and that of ip{u) +7 at the corresponding points u = and u = 1 are the same. 



From (13.151) . we have 

dw dip /df(t (w),w) 

(3.18) 



du du I dw 

i](u — u + )(u — uJ)w{\ — w) [(2a — 1) — Vl + 46 2 w 2 — Ab 2 w 



-2b 2 u(u — l)(w — w + )(w — w_) 



for u 7^ u±, and by l'HSpital's rule, 



(3.19) 



dw 
du 



yjrf + Aar}{\ - 2a)(l - a)(-r)) 



1/2 



b 3 Vb 2 - 4a + 4a 2 

for u = u±. (Note that a/1 + 46 2 u>| - 4b 2 w T = 1 - 2a). 

By coupling (13. 8 p and (13 . 15j) . the double integral representation of t n (x,N + 1) given 
in (12.61) is now reduced to the canonical form 



Ux,N + r 



(3.20) 

where 
(3.21) 



-l) n T(n + A^ + 2) 
2vri r(n + l)r(iV-n + l) 
"+°° ' h(u,r 



3 JV 7 



X 



00 J-y u 



U — 1 



h(u, r) 



u — 1 dto dt 
w — 1 du dr 



and 7„ is the steepest descent path depicted in Figures [6] and [7| when a_ < a < | and 
< a < a_, respectively. 



Before concluding our discussion of Case (i): < a < |, let us give a summary 
of the reduction process from (12. 6p to (13.201) . We first deformed the original surface 
[0, 1] x 7! e C 2 in (12. 6 p into p t x 7^ e C 2 , where p 4 is the steepest descent path in 
the t-plane of the phase function f(t,w) in (12.6 j) . which passes through the saddle point 
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to(w) = £q (w) in (13.41) for each fixed w, and where j w is the steepest descent path of 
f(t (w),w) in the u>-plane, which passes through the saddle points w± given in (13.91) . 
Next, for every fixed w, we introduced a one-to-one analytic map t — >• r, defined by (13.61) . 
which maps p t to (—00, +00), and a map w — > u, defined by (I3.15p . which maps 7^ —> j u , 
where 7^ is the steepest descent path in the u-plane for the canonical integral (I3.20p . In 
the sense of Kaminski [18], we may refer to S TXU : = (—00, 00) x 7^ as a steepest descent 
surface in C 2 . Thus, the 2-dimensional transformation (t,w) — > (r,u) defined by 

bln(l-t) + (l-b)lnt + alaw-ahi(w-l) + b]n[l-(l-t)w) 
(3.22) 2 i / \ 

= — r + a In u — a m(u — 1) + r/u + 7 

is, in fact, given in two stages; see (12. 7p . (13. 6 p and (13.151) . This mapping (t, w) — > (t,u) 
takes p t x 7^ in the C 2 -plane (t, w) into a steepest descent surface S TXU in the C 2 -plane 
(r, u). Hence, in conclusion, we have made the following change of surfaces of integration 

[0, 1] x 7l e C 2 txw — ► Pt x 7w e C 2 XW — > 5 rxu e C 2 XU . 

Case (ii) —00 < a < 0. 

In this case, instead of (12.61) . we use the double integral representation (12.91) . The 
argument here is similar to that in Case (i). But, we note from §2] that when a is neg- 
ative, both W- and w + are real and W- < < 1 < w + . Taking account of the shape 
of the contour 72 in (I2.9p . we need be concerned with only the saddle point W-, and not w + . 

Since a is negative and hence less than ~, the relevant saddle point is again t (w) = 

to (w); see the argument following (13.41) . except with / replaced by /. Thus, (I2.9p can be 
written as 



2ni T(n + 1)T(N -n + 1) 

( 3 - 23 ) r+J r 1 - M 

_J_ e N-f(t ( w ), w ) e -NT*M dwdT 

00 ^72 w - 1 dr 



where 7 2 is the integration path of w in (I2.9P ; cf. (13.81) . However, the mapping w — > u 
differs from that in Case (i). It can be shown that the second derivative of the phase 
function f(t (w), w) is positive at w = u>_, and the steepest descent path passing through 
U7_ is perpendicular to the real axis; see Figure [H 
Recall the Hankel integral for the Gamma function 

1 1 /' (0+) 

(3.24) / e « u -du, 

T(z) 2m J_ 00 

where the contour is a loop starting at —00, encircling the origin in the counterclockwise 
direction, and returning to —00. With z = —aN + 1 and u replaced by — u, we obtain 

1 A[aN f (0+) 1 

(3.25) = — / -e N ^du, 

V ; T(-aiV + l) 2m u 
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Figure 8: The steepest descent path of w, when a < 0. 



where 

(3.26) ip(u) = a ln(— u) — u, 

and the saddle point of if) is at u = a. 

Observe the similarities of the saddle points and the integration paths (see Figure | 
in the two integrals (13. 23j) and (13.251) . Make the transformation w — » u defined by 

= aln(— u) — u + 7 

with 

(3.28) u(wJ) = a, 

where 7 is a constant to be determined. 

Coupling fl3T23|) and (13371) gives 

-1)" T(n + N + 2) 



Ux,N + l) 



2m r(n + l)r(JV - n + 1) 



( - 3 ' 29 ' ) r + ™ /-( 0+ ) h(v t) ~ 



00 J +00 " 



where 

u dw dt 



(3.30) h(u,r) = 

w — 1 du dr 

From (I3.27p . we have 

dw a — u I df(t (w),w) 

du u / dw 

(3.31) / 

(a - u)w{l - w) [(1 - 2a) + y/l + 46 2 w 2 - 4o 2 w 



2ub 2 (w — w + )(w — w_) 

for u a, and by l'Hdpital's rule 



dw f (l-a)(l-2 



a 



< ^ 3 ' 32 ' ) dw \ 6V& 2 - 4a + 4a 2 

for u — a. 



1/2 
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4 THE MAPPING w -» u IN ( 137151 ) 

Before establishing the existence and uniqueness of the constants rj and 7 in equation 
(13. 151) which defines the mapping w — » u, we first need a preliminary result. 

Lemma 1. Let < a < \, and put g(rj) := ip{uJ) — ip(u + ) when a_ < a < |, where 
u± — (rj ± \/r] 2 + Aar])/2r] are given in ^3. 14\ ); put g±(r]) := ip(u_ ) — ip(u+ ) if/ten 
< a < a_ ; where u± denote the u± on the upper edge of the real line and denote 
the u± on the lower edge of the real line; see Figure^ We have 

(i) when rj < —4a, g+{ri) = g~{rj) = k(rj), where k{rj) is a real and monotonically 
increasing function with range (— 00, 0]; 

(ii) when —4a < rj < 0, g(rf) is purely imaginary and —ig{rf) is a monotonically increas- 
ing and continuous function with range (— 2a7r, 0). 

Proof. In case (i), we have rj 2 + 4arj > 0. By definition, 

9±{v) = aln it_ — aln(u ( L' 1 — 1) + rju^ — aln-u^ + aln(u^ — 1) — rju+ . 

Recall the movements of u± shown in Figure [7J Since u^ ± ^ lie on the lower edge of the cut 
[0, 1], it follows that 



«!">-! = I*," - lie- = ^-V^ 2 + 4a V „ 
1 2t] 



Similarly, since lie on the upper edge of the cut, we have 



r4 +) - 1 = - lie" = 7/+v/r?2 + 4a V. 



k(rj) := 2a In- + ~JH_ _ ^ 2 _|_ 4^ 



Put 

fcf-n") ■= 9n In ^ 

rj + \J rf + Aarj 
Then 

= ^-(^) = 



Let £(77) = \/r] 2 + Aarj/rj. It is easily checked that £(r]) is a monotonically increasing 
function in —00 < r] < —4a with range (—1, 0]. If 

Jfe(f) := 2aln^^ + 4a : ; 



then, 



i.e., fc(£) is monotonically increasing in — 1 < £ < 0. Since £(77) = /c(£(r/)), we conclude 
that fc(^) is monotonically increasing in —00 < i] < —4a with range (— 00, 0]. Note that 
k(rj) = when rj = —4a, i.e., a = a_. 
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In case (ii), we have rj 2 + Aar] < 0. Hence, 



u± 



± iyj —Aar] — r] 2 
2~q ' 



From (13.131) . it is clear that Re ip{u + ) = Re ip{uJ) and Im i/j(u^ 
g(rj) is purely imaginary. Since rj is negative, we have 



i21m 



ri — i\Z—ri 2 — Aar] —n — i\/—r] 2 — Aar] 
a In a In 



2r] 



(■ 



^-rf - Aar] 
2a arctan tx + arctan 



2rj 

\f—rj 2 — Aar] 



-Im ip(uJ). Thus, 



— iy—rf — Aar/ 



-v 



—r] 2 — Aar] 



Since —ig'(rj) = —yj—rf' — Aarj/r) > 0, the function —ig{rj) is monotonically increasing. 
Furthermore, since g{— 4a) = —2nai and g(0) = 0, the range of this function is (— 2an, 0). 



□ 



Theorem 1. Let < a < |, (i) When a_ < a < |, there exist unique real numbers rj 
and 7 satisfying the system of equations u{w + ) = U- and u(w-) = u + in 113.16]) . i.e, 



(4.1) 



b ln(l - t+) + (1 - 6) In t+ + a In w + - a \n{w + - 1) + b ln[l - (1 - t + )w^ 
=a In m_ — a ln(u_ — 1) + r]U- + 7, 
(1-6) lnt_ + alnw_ - aln(u>_ - 1) + 61n[l - (1 - t_)u>_ 
=a In u + — a ln(w + — 1) + r]u + + 7; 



61n(l - f_) 

see a/so (fi?.7| ) and Ii3.15\) . (ii) When < a < a_ ; i/iere exzst unique real numbers r] and 7 
satisfying the system of equations in \3.17a\) ( since §3. 1 7b) is equivalent to \3.17a\j ), i. e, 



61n(l - t+) + (1 - 6) lnt + + alnw^ - a\n(w ( +> - 1) + 61n[l - (1 - t+)w 



(+) 



,(+)l 



(4.2) 



-a\nu 



(+) 



aln(«L +) - 1) + V u{ -' + 7 



(+) 



61n(l + (1 -6) lnt_ + aln«/ +) -aln(w w - 1) + 61n[l - (1 -t_)wi" 



(+) 



,(+)l 



=alnw 



(+) 



aln(w+~^ — 1) + r]u^' + 7 



(+) 



Proof. When < a < a_, we denote the left-hand sides and the right-hand sides of the 
two equations in (14.21) by L±, L2 and Ri, R2, respectively. We recall that both w±~^ are 
real and lie in the interval [0, 1], and that they move along the upper edge of the real line; 
see Figure [51 Thus, in L\ and L2, we have 



aln(wi +) - 1) = a\n(l - w ( ±>) + am. 

In this case, since < t± < 1 and < w± < 1, it can easily be verified that 1 — (1 — 
t±)w± > 0. Hence, the terms in L\ and L 2 are well defined. Furthermore, since all terms 
containing t± are real, we have L\ — L2 = Re L\ — Re L 2 . Also, for < w < 1, 

Re f(t (w),w) = Mn(l - t (w)) + (1 - b)lnt (w) + alnw 
-aln(l - w) + 61n[l - (1 - t (w))w], 



(+)^ 



15 



which is a real-valued function. It can be shown that Re f(to(w),w) attains its maximum 
value at w — w_. Hence, Re /(t_,u;_) > Re f(t + ,w + ). and there exists a non-positive 
number r(a) such that L\ — L 2 = r(a). By Lemma []Ji), g+{i]) — where k(r]) is 

a real- valued monotonically increasing function with range (— oo,0]. Therefore, for each 
a G [0, a_], there must exist a unique r\ such that k(rj) = r(a) and the equations in (I4.2p 
hold. The value of 7 is then determined by either one of the two equations in (I4.2p . 

When a_ < a < ~, we also denote the left-hand sides and right-hand sides of the 
two equations in (14. ip by L%, L 2 and Ri, R 2 , respectively. Recall that w± are complex 
conjugate numbers. Since f(t,w) = f(t,w), from (14. ip we have L\ — L 2 = q(a)i, where 

q(a) = 21m {61n(l - t+) + (1 - b) lnt+ + alnw+ 
-a\n(w + - 1) + 61n[l - (1 - t+)w+]}. 

Straightforward calculation gives 
q(a) = 



-y/Aa - 4a 2 - b 2 ^ 6a/4o - 4a 2 - b 2 

2b arctan h (1 — 0) arctan ■ 



2 -2a + 6 v ' 2-2a-b 2 

\/4a - 4a 2 - 6 2 

+2a arctan air 

b 

Note that in the present case, 4a — 4a 2 — b 2 > 0. Also, we have 



, A/4a - 4a 2 - b 2 
q (a) = —2ix + 4 arctan < 0. 

Thus, q(a) < o(a_) = — 27ra_ < 0. Now, consider the function g(a) + 27ra. Since 



d . . . . V^a - 4a 2 - b 2 

— (q(a) + 27ra) = 4 arctan > 0, 

da b 

it follows that q(a) + 27ra > a(a_) + 27ra_ = and hence g(a) > — 27ra. Coupling these 
two results, we obtain — 27ra < q(a) < 0. From the right-hand sides of the two equations 
in ( 14. ip . by Lemma [T](ii) we have R\ — R 2 = where r(r]) is a real-valued increasing 

function with range (— 2a7r, 0). Therefore, for every fix a G (a_, |], there exists a unique 77 
such that r(r)) = q(a) and the system of equations in ( 14. II) holds. The value of 7 is again 
determined by either one of the two equations in (14. ip . □ 

Theorem 2. The mapping w — )■ u defined in Ii3.15\) is one-to-one and analytic when 
< a < \. 

Proof. As in the cases of Charlier polynomials [7] and Meixner polynomials [8], we intro- 
duce an intermediate variable Z defined by 

^ 61n(l - t ) + (l - b) Into + a\nw - a\n(w - 1) + 61n[l - (1 -t )w] 
=Z = a In u — a ln(u — 1) + f]u + 7, 

where t = t (w) is the relavant saddle point t^ (w) given in (13.41) . To establish the prop- 
erties of the mapping w — > u, we again divide our discussion into two cases: < a < a_ 
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Figure 9: The upper half of the w-plane. 



(1 



and a_ < a < 2- 

When < a < a_, we consider the upper half of the tu-plane. Since the functions 
f(t (w),w) and ip(u) + 7 are both symmetric with respect to the real line, the case for 
the lower half of the w-plane can be handled in the same manner. To avoid the multi- 
valuedness around the saddle point, we divide the upper half of the w-plane into two 
parts by using the steepest path through ; see Figure [9j Call these parts region I and 
region II. Note that there are branch cuts, one along the infinite interval (—00, 1] and the 
other along the vertical line Hew = ~ from {b + — b 2 ) /2b to | + ioo. The cut along 
the vertical line was introduced in (13 .4ft . 

As w traverses along the boundary ABB\CC\DD\EFGG\HA of region I, the image 
point Z traverses along the corresponding boundary of a region in the Z-plane; see Figure 
[TOl In Figure [121 we draw the boundary of the region, corresponding to region I, in the 
w-plane. The image point Z = i^(u) + 7 traverses along the boundary of the same region 
shown in Figure [TU1 



A 



IC. 



11 



L-1 



,4 



^ It 1G IL 

Figure 10: The image of Region I in the Z-plane. 
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Figure 11: The image of Region II in the Z-plane. 




A R >, A L K Ci 5/ JjJ J ui Cf 

Figure 12: The upper half of the it-plane. 



Let tp(w) = f(t (w),w) so that ( 14. 3 ft becomes <p(w) = ip(u) +7. Since the mapping 
w — y u is just the composite function ip~ l and ip(u) + 7, it is one-to-one on the boundary 
of region I. By the same argument, one can prove that this mapping is also one-to-one 
on the boundary of region II. For the image of region II in the intermediate Z-plane, 
see Figure [HI By Theorem 1.2.2 of [20l p. 12], the mapping w — > u is one-to-one in the 
interior of both regions. As explained earlier, the one-to-one property of this mapping 
in the lower half of the w-plane can be established by using symmetry of the functions 
with respect to the real axis. Note that the only possible singular points of the mapping 
w — y u are at w = 0, 1 and w±~\ Since the images of these points in the w-plane are 



bounded (see the points 0,1 and in Figure [12]), the mapping is indeed one-to-one and 
analytic in the whole u>-plane. 

When a_ < a < |, the one-to-one and analytic properties of the mapping w —> u can 
be proved in the same manner as in the previous case. Hence, we simply skip this part of 
the argument. 

Thus far, we have proved the one-to-one and analytic nature of the mapping w — > u 
defined in ( 13.151) when < a < a_ and when a_ < a < \. To show that it also has these 
properties at a = a_, we note that equation (I3.15P can be written as 7, /) = 0, 
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where g is just the difference between the left-hand side and the right-hand side of the 
equation. Since rj, 7 and f(to(w),w) are all continuous in a (see the proof of Theorem 
1), by the implicit function theorem u is continuous in a except at u±, where dg/du 
vanishes. The continuity of u in a when u = u± comes from (I3.14p and the continuity of 
i]. Therefore, the mapping is also one-to-one and analytic when a = a_. □ 



5 ANALYTICITY OF h(u, r) 

Before deriving the asymptotic expansions of the integrals in (I3.20p and (I3.29p . we first 
investigate the analyticity of the amplitude function h(u, r) in the neighborhood of the 
steepest descent surface S TXU = (—00, 00) x 7 n ; cf. the statements following (I3.22p . For an 
expression of this function, see (I3.2ip for the case < a < |, and (13.301) for the case a < 0. 

Case 1: < a < ~. 

For w ^ w± and t ^ t+{w) (i.e., r ^ 0), we have from fl3T2T]) . f l3TT8|) and (jBTTj) 
_ u- 1 2ri(l - (1 -t)w] 

i](u — u + )(u — u)w(l — w) [(2a — 1) — \/T+ Ab 2 w 2 — Ab 2 w 
x 



—2b 2 u{u — l)(w — w + )(w — W-) 

Since p t x j w in the t x w-plane does not pass through t$ {w) (unless £q (w) coalesces with 
to (w), a case which we will discuss next), we may substitute (13.41) in (15. ip to give 

2wrt(l - t)[l - (1 - t)w] 

h(u,T) - 



[2(1 + b)wt - 2w + 1 + ^l + Ab 2 (w-l)w]{t - t$(w)) 
r]{u — u + )(u — uJ)w [(2a — 1) — VT+ Ab 2 w 2 — Ab 2 w 



b 2 u(w — w + )(w — W-) 

Thus, in the neighborhood of the steepest descent surface S TXU , h{r,u) is analytic except 
possibly at points u = and u — 1, which correspond to w = and w — 1, respectively. 
From (13 .4p . it is easily verified that £q ( w ) — > 1 — b when to — >• and £q (u;) — > i when 
if — )■ 1. Substituting these values into ( 13 . 1 5|) or (I3.22p yields — = 0(1) as w — )■ and 
0(1) as w — y 1. Hence, neither u = nor it = 1 is the singular point in this case. 



w— 1 

u-1 



If w = w± but t 7^ £q (w), then in (I3.2ip we use (I3.19P instead of (I3.18P to arrive at 

u T - 1 2rt(l - t)[l - (1 - t)w±] 



h(u, t) 



w ± -l(l + b)w±{t - t±){t - t (w ± )) 
vV + 4arj(l - a)(l - 2a)(-T]) ^ 



b 3 Vb 2 - 4a + 4a 2 

Since we are in the case < a < |, we have 7^ 1 (i.e., u + 7^ 1), to + 7^ (i.e., u_ 7^ 0) 

and a 7^ |(1 + yjl — b 2 ). Thus, w + = 0, U- = 1, and a = 1 ~ v 2 1 ~ j ' 2 are the only possible 
singularities, which by (I3.15P and (I3.16P correspond to W- = 0, w + = 1 and rj = —4a, 
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respectively. It can be also verified that ^- = 0(1) as u + — > 0, ^ + _\ = 0(1) as ?i_ — > 1, 

and a/// 2 + Aar]/Vb 2 - Aa + Aa 2 = 0(1) as a — y |(1 — a/1 — b 2 ). Hence, these possible 
singularities are actually removable. 

If w 7^ w± but t = tg (w), then in (I3.2ip we use the second equation (instead of the 
first one) in (13. 7p to reach 

2t +H(i-t + H)[i-(i-t+HH ' 1/2 

a/1 + Ab 2 {w - l)w 
i](u — u + )(u — u_)w \{2a — 1) — a/I + 46 2 w 2 — Ab 2 w 

x 



/i(w, r) 



2b 2 u(w — w + )(w — u>_) 

where, in addition to the possible singular point at u = which we have just discussed, 
we have introduced two new possible singularities w = (b ± i\Jl — 5 2 ) /26. However, from 
Figured] it is clear that these two points will not be in the neighborhood of ptX'jw, unless 
a — > ~ in which case w — > w± (see the discussion below). 

If w = w± and t = ^(w), then in (" 13 . 2 1 [) we use (13 . 19[) and the second equation in (13.71) . 
Since £q (u>±) = t±, the result is 

/ 2t±(l - t±)[l - (1 - 1±)«;±] 

n{u,T) = <^ 

w± — 1 [ b d 

(5.2) . x 1/2 

a/t? 2 + Aar](l - a)(l - 2a)(-r?) 1 

a/& 2 - 4a + 4a 2 

Note that a = ~ is not a singular point, and since \/r] 2 + Aarj/y/b 2 — 4a + 4a 2 = 0(1) 
as a — > |(1 — a/1 — b 2 ) as we have discussed before, there is no singular point in the 
neighborhood of the integration surface. 

Case 2: a < 0. 

In this case, only the negative saddle point u>_ is relevant, and the steepest descent 
path associated with w_ will not pass through w + . Thus, the situation is much simpler 
than the case < a < ~. By using (I3.30p . (13. 7p . (13 . 3 1 [) and (I3.32p . we have the following 
analytic expressions for h(u,r): 

If w 7^ and t ^ ^(w), we have 

(u - a) [{I - 2a) + a/1 + Ab 2 w 2 - Ab 2 w\ rt(l - t) [1 - (1 - t)w\ 



h(u, t) 



b 2 (w - w+)(w - wJ) (1 + b)(t - U(w))(t - t (w)) ' 

Since w ^ W- and t ^ t$ (w), clearly there is no singularity in the neighborhood of the 
integration surface. 

If w u>_ but if: = (w), then 
h(u, r) = 



2t+{w){l - t+{w))[l - (1 - 4{w))w] 1 V2 (u - a)w [(1 - 2a) + a/1 + Ab 2 w 2 - Ab 2 w 



a/1 + 46 2 (w - l)w j 2b 2 (w -w+)(w -w-) 

20 



Again, there is no singularity. 

If w = u>_ (in which case u = a by (I3.28P ) but t ^ £q (w), we have 

-2rt{l -t)[l - (1 -t)w-] f ab(l-2a) 1 1/2 
~ (1 + 6)(t-t_)(t-to(w_)) 1(1 -a)Vb 2 -4a^4^ J 

The only possible singularity occurs when if_ = (i.e., when a — > 0~). Since it can be 
shown that — = 0(1) as a — > 0~, h(u, r) is again bounded and analytic as a — > 0~. 

If w = u>_ and t = td(w), then 
(5.3) _ J " *-)[! " (1 " *->--] 1 V2 



(1 - w-)by/b 2 - 4a - 4a 2 
and there is no possible singularity at all. 

The conclusions in all these cases infer that for a < 0, h(u, r) is analytic in the neigh- 
borhood of the surface of integration. 

6 ASYMPTOTIC EXPANSIONS 

Case 1: < a < \. 

Let ho(u,r) = h(u,r), where h(u,r) is the amplitude function given in (I3.2ip . Write 
(6.1) h (u, t) = a (r) + b (r)u - (u - u + ){u - u-)g (u, r). 



Clearly, 



(6.2) 



and 



u-h (u + ,r) - u + h (u-,r) 
ao{r) = 



bo(r) 



u_ — u + 
h (u_,r) - h (u + ,r) 



U- — «j 



(6.3) g (u, t) = i [h (u, r) - a (r) - b (r)u] . 

{U — U + ){U — U-) 

For u 7^ u±, go(u, r) is analytic in u wherever ho(u, r) is analytic. Furthermore, it is easily 
shown that the limits 

, , h' Q {u + ,r) h (u + ,r) - h (u^,r) 
hm g (u,T) = - if a^a_, 

u-*u+ U + — U- \U + — U-Y 

h' (u- } r) h (u + ,r) - h (u^,r) 
hm g {u,T) = 1 - if a^a_, 

lim g (u,T) = lim. g Q (u,T) = -hl(u-,r) if a = a_, 
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exist and are finite. Thus, go(u, r) is actually analytic in u everywhere ho(u, r) is analytic. 



Recall the integral representation of the confluent hypergeometric function M(aiV + 



(6.4) 



M(aiV + 1, 1, r/N) = -L I e 



dw 



where ip(u) = ahxu — aln(u — 1) + rju and 71 is the loop contour which starts at u — 0, 
encircles u — 1 in the positive direction, and returns to u = 0; see (I2.5p . From (13. lip , we 
also have 



(6.5) 



M'(aN + 1,1, r]N) 



1 
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27U ./ 7l u — 1 



where the derivative is taken with respect to the third variable in the function M(c, d, z). 
Inserting (16. ip in (I3.20p . we obtain 



(6.6) 

t n (x,N + l) 



(-l) n T(n + N + 2) 
T(n + 1)T(N -n + 1) 



D N 1 



a (r)e- NT2 dT 



M(aN +1,1, r]N) 

/+00 
6 (r)e- 7Vr2 dr + e+ 
■00 



—00 
+00 



where 



: i 



2ni 



(u - u+)(u - u-)g (u, r)e 



du 



e~ Nr dr. 



To the inner integral, we apply an integration-by-parts; the result is that it is equal to 

, u=0+ 



r]N 



ue N ^g {u,r) 



u=0- 



1 



/ il ( M ,r)e^ w 



du 

u-r 



where 



hi(u,r) = [u - 1) 



g {u,T) +u 



dgo(u,r) 
du 



By using (I6.3p and the fact that ho(u) = 0(1) as u — > as we have shown in §3 the first 
term above is actually zero. Thus, 



(6.7) 



2-Kir]N 



OO J~j u 



u-1 



e- Nr2 dudr. 



The two integrals in (I6.6P can be asymptotically evaluated by using Watson's lemma 
p. 20]. So, let us expand a (r) and b (r) into Maclaurin expansions 



j=0 



j=0 



r+ 
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Since 
(6.8) 

we obtain 

t n (x,N + l) = 



/■ +00 r 2 Je -^ dr= r (j + | 

J — OO 



/+oo 
T 2 j+ l e -Nr* dT = Q) 
-oo 



■i)-r(» + j V + 2) .^, M(ajV + Ui)M 



Co 



N 
d 



+ e 



T(n + 1)T(N -n + 1) 

+M' (aN + 1,1, rjN) 

where Sf and Sf denote error terms that are 0(N~ 3 / 2 ). 

The above procedure can be repeated. For I — 0, 1, 2, • • •, we define recursively 
(6.9) 
(6.10) 



i r ' 



hi(u,r) = ai(T) + bi(T)u-(u-u^)(u-u + )gi(u,r), 

dgi(u,r) 



h t+1 (u,T) = (u-1) 



gi(u,r) + u- 



du 



and 
(6.11) 



j=0 j=0 



Using induction, one can show that hi(u,r) and gi(u,r) are both analytic and bounded 
as u — y 0. Hence, 



(6.12) 
t n (x,N + l) 



{-l) n T{n + N + 2)e N ^ f 



T(n+1)T(N -n + 1) 1 



M(aN + 1,1, rjN) 
+M'(aN + 1,1, rjN) 



.1=0 iv 



p-1 



z=o 



+ e 



v f ' 



where 
(6.13) 

(6.14) 
(6.15) 



\^ 2m + 1 

Q — 2_^ ai - m ' 2m v — ^ — J' 



m=0 



— ^ &I-m,2 m r( ), 

m=0 



,+ 



— — / / h p (u,T) -e T dwdr, 

irjNNPj^ L IA u-1 



2ixinN Np J_ 



oo J 7„ 



and 5+ 5+ are error terms of order 0(N p 1 / 2 ). 
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To prove that for fixed b = n/N e (0, 1), (I6.12p is a uniform asymptotic expansion in 
a G [0, |], we need to show that there are positive constants A p and B p , independent of 
a, such that 

A B 
(6.16) |e+| < — ^r|M(aJV + l,l,77iV)| + — ^ I M' (aN + 1, 1, nN)\. 

p N p+ ^ 



Case 2: a < 0. 



Proceeding in the same manner as in Case 1, we put ho(u, r) = h(u, r), where h(u, r) 
is given by (I3.30p . and write 



(6.17) 



h (u,r) = a (r) + (u - a)g (u,T). 



Since lim g Q (u,r) = h' (a,r), g (u,T) is analytic in u wherever h (u,T) is. By (I3.29P and 

u— ¥a 

(I3.25p . we have 



(6.18) t n (x,N + l) 



where 



aN N-y 



(-l) n T(n + N + 2)N~ aN e 
T(n + 1)T(N -n + l)T(-aN + 1) 



+oo 



ao(r)e-^dr + er 



N an T{-aN + 1) 
2th 



+oo 



(0+) _ J..' 

(u-a)g (u,r)e N ^ d " 



+oo 



U 



e-^ 2 dr 



and V , ( M ) is given in (13. 26 p . To the inner integral, we apply an integration- by-parts. The 
result is that it is equal to 



Nip(u) 



+ 



N 



/i 1 (u,r)e iV ^ ) 



+oo 



du 



u 



where 



hi(u, t) = u 



dg (u,T) 
du 



Since ip(u) ~ — u as u — > +oo, the first term vanishes and we have 
N aN T(-aN + 1) r+ °° 



2niN 



(o+) _ d . 

/n(ii,T)e*W«> — 



dr. 



The Laplace integral in (16. 18ft can again be asymptotically evaluated by using Watson's 
lemma, and the result is 



(6.19) t n (x,N + l) 



aN p N-f 



(-l) n T(n + N + 2)N~ aN e 
T(n + 1)T(N -n + l)T(-aN + 1) 



^0,0 r A\ ,-,<:- 



where ao,o is the constant term in the expansion of Oo(t) and 5 1 denotes the error term 
resulting from Watson's lemma; see (I6.8p . 
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The process can be repeated. For I = 0, 1, 2, • • we define 
(6.20a) 
(6.20b) 



h t (u, t) = ai(r) + (u- a)g t (u, r), 
dgi(u,r) 



h l+1 (u,r) = u- 



du 



+oo 



(6.20c) 

As before, we arrive at 



j=0 



(6.21) 

where 

(6.22) 

and 

(6.23) 



aN 



Ux,N+l) 



(-l) n T(n + N + 2)N~ aN e 
T(n + 1)T(N -n + l)r(-aiV + 1) 



\^ 2m + 1 

Q — 2_^i a l-m,2rnL { ^ ) 

m=0 



p-1 



1=0 



\- Ci 



N aiy T{-aN + I) 
2niNP 



+oo 



(o+) „ ch : 

h p (u,r)e N ^ — 

+oo u 



e p +S p 



e-^ 2 dr. 



The error term 5 p comes from applying Watson's lemma to the functions a;(r), < / < 
p — 1, and is of order 0(N~ P ~^). Thus, to prove that expansion (16.211) is uniform with 
respect to a, we only need to show that there exists a constant M p , independent of a, 
such that 



(6.24) 



Ie~| < 



M p 



7 ERROR ESTIMATES 

To estimate the error terms e p and e~, we shall actually provide two slightly stronger ver- 
sions of the results stated in (16. 16ft and (16.241) . Indeed, we shall show that the constants 
A p and B p in (16.161) are independent of b 6 [5, 1 — 5} and a e [0, |] for any 5 > 0, and 
similarly that the constant M p in (16.241) is independent of b G [5, 1 — 5} and a G (— oo, 0). 
That is, the expansion in (16.121) is uniform in b G [5, 1 — 8] and a G [0, |], and the expansion 
in (I6.2ip is uniform in b G [5, 1 — 5} and a G (— oo, 0). 



To proceed, we will divide our discussion into five cases: 

(i) a_ < a < \ and N% (rj + 4a) is unbounded; 

2 

(ii) Ns^r] + 4a) is bounded; 

(iii) < a < a_ , N 3 (j] + 4a) is unbounded and aN is unbounded; 

(iv) a > and aN is bounded; 
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(v) a < 0. 

From (I3.15p . we note that a = a_ corresponds 77 = —4a. 

Case (i) First we recall the following result of Fedoryuk [17, p. 124-125]; consider the 
integral 



/(A) 



f(z)e xs{z) dz, z = ( 



Zl, 



where A > 0, dz = dzi ■ ■ ■ dz n and 7 is a smooth manifold of (real) dimension n. Assume 
that the function f(z) and S(z) are analytic in some domain D containing 7. A point z° 
is called a saddle point of S(z) if VS(z°) = 0, and a simple saddle point if det S"(z°) 7^ 0. 

Theorem 3. Let max i?e S(z) be attained only at a point z° , which is a simple saddle 
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point and an interior point of '7. Then, as A — > +00 ; we have 



2n 



/(A) = ( T ) [det S 



n/2 



/// On] -1/2 AS( 2 °) 



/(: n ) + 0(i 



T/ie choice of the branch for the root depends on the orientation of the integration surface. 

In our present case, n = 2 and there are two simple saddle points (u + , 0) and 0). 
By Theorem [31 it follows from (I6.15P that as iV — > 00 



6 ~ 
P 



(rt + Aa)- 1 /^/ 2 
V2i(-7]) 5 / 4 NP+ 2 



h p (u + , Q^ e i ( aN ^- 2aN 'P-+^ N V- r i 2 - 4ar i + ^—^-) 
+h p (u_, {^ e - i i aN ^~ 2aNl P-+l N y/- r i 2 - Aar i + '^f-v-) 



where ip_ is the argument of w_ in (— tt, 7r]. The asymptotic formulas of lSA(aN + 1, 1, 77 
and M'(aAf + 1, 1,7/iV) can be obtained in the same manner with n = 1 and h p (u,r) 
replaced by 1. Thus, 



M(aN + 1,1, r)N) 



(r/ + 4a)" 1 / 4 e iV ^ 2 



v / 2TiVi(-?7) 1 /4 



3 i(aAr7r-2aA r <p_ + |Af A /- J? 2 -4ar7+^-ip_) 



+e 



-i(aNTr-2aN<p- + ±N^-r) 2 -4:ari+22—ip„) 



i(aNTT-2aN<fi- + ±N^-r] 2 -4:ari+^--2ip-) 



_^_ e -i{aNw-2aNip^ + ±Ny/-T 1 2 -4:ar ] +?l-2 l fi-) 

A combination of the last three equations gives the error estimate in (I6.16P . 

Case (ii) Here, N 2 ^ 3 (i] + 4a) is bounded. Hence, 77 — > —4a and, correspondingly, 
a — > a_; see the statement preceding (13.151) . For a < a_, the contour 7 n in (16.151) is shown 
in Figure [3 for a > a_, the contour 7^ is depicted in Figure [6j We divide 7^ into two 
parts, so that the error term e+ in (16.151) can be written correspondingly as e+ = + e+ 2 . 

We call the part of 7„ in the upper half of the plane 7^ , and the part of 7„ in the lower 
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half of the plane 7„ . Because of the cut in the interval (0, 1), the phase function ip{u) of 
the u-integral in (I6.15P can be written as 

ip(u) = a In w — a ln(l — u) + rju — ian = ip{u) — ian 

for u G 7« 1 , and 



ip(u) = ip(u) + iair 



for u £ 7u . Since ^(w) has two coalescing saddle points w-t given in ( I3.14p . we make the 
standard transformation 

(7.1) ^ U) = l v z-Cv + A 

with the correspondence of the saddle points u± and ±a/C given by 

(7.2) M_ ^ y/£, U + ^-^. 

If C > 0, then are both real; if ( < 0, then i^/C are complex conjugates and purely 

imaginary. The values of A and ( can be obtained by using (17. ip and (17.21) . The images 
of the contour 7 ; |, 2 ' > under the mapping u-^t) defined in ( 17. ip are depicted in Figures 
and [HI The part of the error term e+ now becomes 



g— i<nrN+^r)N poo 



e p,l 



2mr ] NP+ 1 J^Jc, 
where 

1 du 



(v,r)e m ^-^e- NT2 dvdr, 



(7.3) g(v,r) = h p (u,r) 



u — 1 dv 

















> 



Is/ 



Figure 13: Contour C 2 (a < a_). 



Figure 14: Contour Co (a > a_). 



The part 2 can be handled in exactly the same manner, and we have 

g(v, T)e N ^ vi -^e- Nr2 dvdT, 



^ianN+^rjN poo 



~ p ' 2 2mr]NP +1 



oo JC 3 



where the contour C3 is the reflection of C 2 in the f-plane with respect to the real axis. 
Since the two contours C 2 and C3 together are equivalent to a contour running from 
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ooe l7r / 3 to ooe J7r / 3 in the right half-plane, in terms of Airy integral Ai(-) we have by the 
Bleistein method J2TJ p.366-370] 



(7.4) 



e; ~Re 



( -iairN+±r)N 



-goe^W-iAi' f N^e 



where 

(7.5) 

(7.6) 



^(v / C>o)-( 7 (-v / C,o) 



For a detailed discussion of uniform asymptotic evaluation of integrals, see [21] Chap. 7]. 

Following the same argument, one can derive the asymptotic approximations of the 
confluent hypergeometric functions M(aiV + 1,1, rjN) and M'(aiV +1,1, rjN). Indeed, we 
have 



M(aN + 1,1, r)N) ~Re 



(7.7) 



-iairN+±riN 



-gJe^JV-f Ai' f A§C^ 



where £ is again determined by (17. ip and (17. 2p . and is the same as that appearing in ( I7.4p . 
The coefficients and are similar to those given in (17. 5p and ( 17. 6p . except that the 
function h p (u, r) in (I7.3P is replaced by 1. The asymptotic formula of M'(aA + 1, 1, r)N) 
is the same as (17. 7p . only with h p (u, r) in (17. 3p now replaced by u. 



A comparison of (17.41) with (17. 7p and a corresponding formula for ~M.'(aN + 1, l,r]N) 
again establishes ( I6.16p . 



Case (iii) In the present case, the contour 7^ in ( I6.15P is shown in Figure [7J As in 
case (ii), we divide 7„ into two parts, and label the parts above and below the real axis 
by ju and \ respectively. Since 7« and 7! are symmetric with repect to the real 
axis and h p (u,r) = h p (u, r), the error term e+ can be written as 



irr)NP +1 



Im 



00 j 7. 



h p (u,r) 



(2) u — 1 



In view of the equation preceding (17.11) , we can further write 



(7.8) 



— < cos(a7rA)Im 

nr)NP +1 \ v ' 

+ sin (got A) Re 



00 

CO 



f h^^ e NHu) e -Nr^ dud \ 
(2) U — 1 J 

h^ e ^u) e -Nr* dud \ 
(2) 11 — 1 / 



OO J J,, 
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where ip(u) = alnu — aln(l — u) + r\u. 

2 

Now observe that since Ns(^r] + 4a) is unbounded, the rate in which r\ — > —4a is 
controlled; that is, the two saddle point u± of the phase function ip{u) do not coalesce 
sufficiently fast; cf. (I3.14p . As a consequence, it is known that the classical method of 
steepest descent continues to apply. (Note that Nip{u) = N(a\nu — aln(l — u) +r]u). So, 
to apply this method, aN must approach infinity; i.e., x is large.) 

In (17. 8p . there are two terms, and both terms contain the same double integral. Since 
^>(u) is real for u £ (0, 1), to get the real part of the asymptotic approximation of the 

(2) 

double integral we need be concerned with only the horizontal part of the contour 7« . 
Furthermore, since ip{u + ) > ip(u-), the real part of the double integral in ( 17. 8p is asymp- 
totically equal to 



(7.9) pv ^ ' e 



h P (u+,0) _j^ (u+) / 2 7T 



-i/)"(u+) N 



(2) 

Since the integrand on the horizontal part of the contour 7« is real, only the curved part 
of contour 7! contributes to the imaginary part of the double integral. Since the cut for 
ip(u) is now from 1 to +00, we can extend the contour into the upper half-plane by 
its mirror image, and form the integration path V which passes through w_ as an interior 
point. The imaginary part of the double integral in (17. 8p is asymptotically equal to 



Inserting the contributions given in (17. 9p and (17.101) into equation (I7.8P gives 



(7.11) k V 

cos(aiiN)h p (u_,0) N j (u ) I 2 f ( N -i\\ 

2% u--i y -$»{uj) L 1 n 

Again, since ip(u + ) > ip(u_), (17. lip implies 



(7.12) 



+ an exponentially small term 



Asymptotic approximations for M(aiV + 1, 1, rjN) and M'(aiV + 1, 1, rjN) can be ob- 
tained in a similar manner, and we have 



(7.13) 



MiaN + !, l,ijj\r) = 2 tl + 0(JV -: 



+ an exponentially small term 
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e N<p(u+) 

M'(aN + l,l,r)N) = —j^l sm(anN) 



(7.14) 



+ an exponentially small term 



[1 + 0{N- 1 



The estimate in fl6U6|) now follows from fl7TT2|) . ( ESJ) and f[7TT4j) . 



Remark. To 6e more precise, the exponentially small term in (7.13) and Ij7.14 ) ^ s 

Q { e N(i>{u-)-i>{u + ) 



Case (iv) Here, we proceed as in case (iii), and reach the error representation given 
in (17.81) . There are two terms in this equation; we begin with the one involving 



(7.15) 



Go = Re 



As explained in the previous case, for this term we need be concerned with only the 
horizontal part of the contour 7^ . Since aN is bounded, a must tend to zero and r] is 
negative; see the correspondence between a and r\ given in ^J3j Thus, the saddle point u± 
behave like 



2 2V V 



11/ 4a 
2 + 2V 1 + 7 



-- + 0(a 2 ) ^0, 
V 

1 + - + 0(a 2 ) 1. 
V 



As u — > 0, ~ alnu + (a + //)«. This suggests that we make the transformation 

ip(u) = a\nv + (a + rj)v = p(v) 

with u = u + corresponding to v — v — —-^r- Recall that r] is negative and a — > 0; thus, 
v > and Vo — > 0. This says that our saddle point vo coalesces with an endpoint. We 
also note that Nip{u) = Nalnv + N(a + i])v, and that aN is bounded. All these point 
out that the steepest descent method does not apply. Therefore, we rewrite ( 17.151) as 



(7.16) 

where 
(7.17) 
Since 



G 2 ~ Re 



00 p+oc 



00 JO 



l(v, r)v aN e N{a+v)v e- NT2 dvdr 



du 
dv 



g(v,r) 



P"(v ) 



h p (u, r) du 
u — 1 dv 



! u 2 (u - iy 



v=Vo v ij)"{u 
as a — > 0, we obtain by Watson's lemma 

^(0,0) 



w 2 (l -2u) 



->• 1 



(7.18) 



Go 



-N{a + r])} 



aN- 



T ^T(aN + l) [1 + OiN- 1 )]. 
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Next, we consider the first term in (17.81) 
(7.19) G 1 = Im ^J°° J M^ll e *) e -^ 2 dwdr 

As in case (iii), since the integrand on the horizontal part of the contour 7« is real, by 
symmetry we have 



(7.20) G 1 = -( r I ^^e N ^e- N ^dudr 



2i \J-oo JL M - 1 
(2) 

where L is the curved part of 7u plus its minor image in the upper half-plane. Recall 
that there is a cut for ip{u), which extends from u = 1 to u = +oo along the positive real 

line. Since u_ = l + -H ~lasa->l, this saddle point coalesces with the branch 

point at u = 1. For -u near 1, we write u = 1 — v with v being small. Then, 

ijj(u) ~ — alnt> + 7] — (77 + a)f . 

This suggests the transformation 

(7.21) = — a In 2; + 77 — (rj + a)<z = 

The double integral in (I7.20j) now becomes 



(7.22) Gi = — 



if r +, M^)* e «*», e -^ d , dT . 



2« 7_oo 7-oo u - 1 dz 



Since Np(z) = Nalnz — N(r) + a)z + Nr) and iVa is bounded, the method of steepest 
descent does not apply. So we rewrite (I7.22p in the form 

G x = -^e Nr > / / z-^giz^Tje-^+^e-^dzdr, 
^ J -00 J -00 

where 

#0, r ) = Tl~ h p( u ' r )- 

u — 1 dz 

Recall that N(i] + a) is negative. To the inner integral, we can use Watson's Lemma for 
loop integrals [2H p. 20]. Since 



= — 1 and -r^ 

*=o,«=i d2: 



1, 



2 =0,n=l 



It — 1 

we have g(0, 0) = /t p (l, 0) and 

(7.23) G, = -VI, 0) ^ '^Vf ^ [1 + O(iV-)] . 

Coupling ([73HD and (17251) . we have from ([H]) 

iV7r[-iV(a + r/)]^+i [ 1 + °( iV 



-sm(aN7r)h p (0,0)J n T ( aN + V < 
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which in turn gives 

1 f sin(aA%)/i p (0,0) fn 



^T(aN + l) [1 + 0(N-')} 



{70] . " p v Np+i{ 7r[-N(a + V )]« N +^ N 

+ 0{e Nr >) 

Note that 77 is negative. 

The asymptotics of M(aiV + 1,1, r)N) and M'(aiV + 1,1, f]N) can again be obtained in 
a similar manner. We have 

(7.25) M(aN + ^T^ r(aiV + 1} [1 + ° {N ~ 1)] 

+ 0(e N "), 

, M'(aN + 1, 1, V N) = r^^i^i r^ + !)[! + O^ 1 )] 

(7.26) 7rr7[-iV(a + r])] aN+1 K ' l y Jj 

+ 0(e Nr >). 

The estimate flBTTBj) follows immediately from (I723|) . ( TT251) and (OH). 

Case (v) We first consider the subcase in which a is bounded; that is, a 6 (— M, 0), 
where M > is a constant. Since in this case only the saddle point u = a matters (see 
(I3.26P ). by applying the steepest descent method directly, we obtain 

P-27) £; = M^2) ^[l + O^-')]. 

Since h p (u,T) is analytic along the path of integration, for bounded a we can easily find 
a constant M p satisfying (I6.24p . 

Next, we consider the subcase: a — > —00. Motivated by a method of Olde Daalhuis 
and Temme [22] . we introduce the rational functions 

(7.28) Ro(u,w,a) = — !— 

u — w 

(7.29) R n+ i{u,w } a) = — (uR n (u,w,a)) , n = 0, 1, 2, • • •, 

where u,w ^ C, u ^ w and w 2 7^ a 2 . Using an induction argument on n, it can be readily 
verified that there are constants Cy, independent of u, w and a, such that 

n-ln-l-i C--U n ^ 1 ^ j 

(7.30) R n {u,w,a) = wV V 7 wTTZTZT? rrrj, n = l,2, 

^— ' (m — w) 1 Hu — a) n+l 

To extend the validity of the asymptotic expansion (16.211) from bounded a to unbounded 
a, we first define 

(7.31) Po( a ) = min{|u — a\ : u is a singularity of /i ( u )}- 



32 



Since we shall be interested in the property of h p (u, r) in u, we may suppress the depen- 
dence on r and simply write h p (u) for h p (u, r). To obtain an estimate on po(a), we denote 
by u±k the singular points in the w-plane, which are the images of w^e ±2km under the 
mapping (I3.27p - fl3.28p . where w_ is the saddle point given in (13. 9p . Thus, we have 

(7.32) f(to(w-),W-) = aln(— a) — a + 7, 

(7.33) f(t (w.e ±2km ),w.e ±2km ) = a ln(-u ±k ) - u ±k + 7. 

Subtracting (I7.32p from (I7.33P gives 

(7.34) ± 2kir(a + b — l)i = a\ia(—u±k) — u±k — aln(— a) + a. 

For large values of a, numerical calculation shows that the nearest singularities should be 
u±\ which are approximately equal to (3.089 ± 7.461z)a. Hence po(a) ~ —7.75a. (Recall 
that a is negative). [Put u±k = av±k in (17.341) . and have 

, fa + b- 1\ . . 

±2fc7r I 1 i = \n(v± k ) - v± k + 1. 

For large a, the left-hand side is approximately equal to ±2kiri. Numerical computation 
of the resulting equation gives v±i = 3.089 ± 7.4611] 



To estimate the error term e p in (I6.23p . we split the loop contour into two parts C 



and £ 2 , and write 



(735) ^ . « != ^£. i to de ^^ i 

NTj-aN + l) f +oa f ^r) A) ^ 



(7.36) e' 2 = V ' ' / / ■^ll± e ^W e -^ dud 
P ' 2 27TtNP J C2 u 

where C\ = {u G loop contour : \u — a\ < c\a\ 9 }, with < c < 1 and | < 9 < 1. A 
simple estimation shows that e~ 2 is exponentially small in comparison with e~ A ] see [22J 
p. 318-320]. 

Let be a closed contour embracing the curve Ci such that as a — > —00, 

(7.37) length Ti = 0(|a| e ) and distance (Ti,Ci) ~ ci|a| . 

Let f2i denote the domain bounded by Fx, and its closure by Similarly, let T 2 be a 
closed curve embracing C\ and lying inside such that as a — > —00, 

(7.38) length T 2 = 0(\af) and distance (r a , A) ~ c 2 |a| e 

with c 2 < c\. Let fi 2 denote the domain bounded by T 2 , and its closure by f2 2 - Define 

(7.39) h p = sup \h p (w)\. 

It is easily verified that as a — )• — 00 

(7.40) sup \ho(w)\ < co\h {a)\, 
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where cq is another constant independent of a. From (I7.30p . we also have 
(7.41) sup_ \Rp(u,w,a)\<C p \a\-W-V p - e , 

where C p is a constant independent of a. By (I7.40p . ( 17.41 j) and (I7.37p . it follows that 



(7.42) 



2m 



R p (u, w, a)h (u)du 



< C'\a\- {2e - 1)P Mc 



where C' p is also a constant. By Cauchy's integral formula and (I6.20b|) and (I7.29p . we 
have upon integration by parts 



h v {w) = / Rq(u : w, a)h p (u)du 

2m J T 



gp-i(u)d(uRo(u,w,a)) 



(7.43) 



2m 



2m 



[u — a)Ri(u, w, a)g p _i{u)du 



— / Ri(u,w,a)hp-i(u)du- — / a p -xRx 
2m J Tl 2m J Tl 



[u, w, a)du, 



where we have used (I6.20al) to get ft ?) _ 1 (a, r) = a p „i{r). On account of f)7.4ip . the last 
integral is /i p _i( i u)0(|a|~( 2e '~ 1 )). Repeated application of (17.431) gives 

(7.44) h p (w) = ^- [ R p (u,w,a)h (u)du + h p ^O(\a\-^ e -^) + ■ ■ ■ + h O(\a\- {2e - 1)p ) 



2m' j r 

By induction, we have from (I7.44p and (I7.42p 



(7.45) 



hp < C p \a\~ {2e - l)p \h { 



where C p is a constant independent of a. Substituting (I7.45P in (I7.35p . we obtain 

A p \a\- {2e - 1)p \ho(a)\V^ 



(7.46) 



1^1 < ? " 



NP+ 1 / 2 



where A p is a constant not depending on a. From (I5.3p . one can show that ho(a) = 
h (a, t) ~ \/^a as a — > — oo. Thus, to have the right-hand side of (I7.46P independent of 
a, we must choose 6 so that {26 — 1) > i > \ for all p > 1; that is, 9 > |. With this 
choice, we can always find a constant M p , independent of a, for each p such that (I6.24p 
holds. 



8 SUMMARY OF RESULTS 

For a < 0, we have from (I6.2ip and (I6.24p 

(-l) n T(n + N + 2)N~ aN e N ^ 



t n (x,N + l) ~ 



r(n + l)r(iV - n + l)r(-aiV + 1) ^ N l+L 2 



ci 
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as N — > oo, where the coefficients q are given in (I6.22p and the constant 7 is determined 
by (I3T271) and (13T28D . 



For < a < |, we have from ( I6.12p and ( 16.1 6p 

oc 

M(a^ + l,l,^E]^i 



2' 

/ \ -1-1 1 . nT . r\ \ at™ r OO 

q 



. Ar ^ (-l)T(n + + 

1 m + 1)1 (TV — n + 1) 
(8.2) 



+M / (aiV+ 1, 1,T]N)J2 





di 



1=0 



as AT 00, where M is the confluent hypergeometric function defined in (16 ,4p . and M' 
is its derivative taken with respect to the third variable. (Note that the second variable 
in M is a constant.) The coefficients q and d\ in (18. 2p are given in (I6.13P and (|6.14j) . 
respectively. The constants rj and 7 are determined by (I3.15P and (I3.16p . 

For a > |, the asymptotic expansion of t n (x,N + 1) can be obtained by using the 
symmetry relation (11. 8p . 

(8.3) t n (aN, N + l) = (-l) n t n ((l - a)N, N + l). 

The regions of validity of the expansions (18. ip and (I8.2p are, respectively, a < and 
< a < |. However, as we shall show in $91 they can be extended to two overlapping 
regions both containing a neighborhood of a = 0. 

The integral given in (I3.12p for the confluent hypergeometric function is a typical ex- 
ample of Airy-type expansion [STJ Chap. VII, §4]. In fact, we have already illustrated 
this point in ( 17. 7p . and the region of validity of this formula is0<5<a<~. Thus, for 
$ < a < |j one can derive from (I8.2p an asymptotic approximation for t n (x,N + 1) in 
terms of the Airy integral Ai(-) and its derivative Ai'(-). 

We now derive from (18. ip and (18.21) two simple formulas for t n (x,N + 1) when x is 
fixed, which in turn means a = 0(1/N). First, we consider the case when x is negative, 
i.e., a < 0. To calculate the leading coefficient Co in ( 18.11) . we use ( 16.221) . ( I6.20cl) and 
( I6.20al) . Thus, we have Co = cto.or (|) and a ,o = «o(0) = ho(a,0). From ( 15. 3p . it can be 
readily verified that 



h(a, 0) = - v ^j- h \ + 0(a 2 ) zsa^O-. 
Hence, 

/ x a/2(1 - b)Tl , 9x 

(8.4) c = - v b3/2 a + 0(a 2 ) as a 0". 

To find the value of the constant 7 in ( 18. ip . we recall (I3.27P and (13. 28j) . and obtain 

7 = f(to(w-), W-) — a ln(— a) + a. 
Since to(w-) = t_, from (I2.10p and (13.21) it follows that 

(8.5) 7 = 61n6+(l-6)ln(l-6)-2aln6 + 0(a 2 ) as a -> 0". 
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Inserting (JH3J and f 18 . 5 j) in (18. ip gives 
(8.6) t n {x,N + l) = 

for fixed x < 0. 



l) n+1 r(n + AT + 2)A^~ 2:r - 2 

r(iv + i)r(-x) 



1 + 0| 7V 



Next, we consider the case when x is nonnegative; this means that < a < \ and aN 
is bounded. Substituting (17351) and rT7726|) in flO]) gives 

-l) n T(n + N + 2)e N "i f T(x + 1) sinvrx 



f n (z,JV + l) 



(8.7) 



r(n+l)r(JV-n + l) I ViV [-iV(o + r/)] aiV+1 tt 



—Co H — rf 



1 + 



+ 0(e 



The coefficients Co and do can be calculated by using (16.130 . ( 16. 141) . ( 13. 14[) and ( 16. 2ft ; the 
result is 



a , 

-c + -a 
V 



a o,o Vo 

V ' . 

■y/nh (u + ,0). 



From (15. 2p . it can be shown that 

(8.8) h (u + ,0) = - 

Hence, 

(8.9) 

By (13.11) and (13.21) . we also have 



W 2 



+ 0(a) 



as a — y . 



a 2y / 7r?7v / l — b , 
Co + -do = h Oia) as a — y . 

7] b 6IA 



5.10) 



w + = l-- + 0(a 2 ), 



b 2 



+ 0(a 2 



1 



1 



-a + 0(a 2 



t_ = (1 — 6) + (1 -6)a + 0(a 2 ). 



1+6 1+6 

Inserting (I8.10p into (j4.2j) . and subtracting the two resulting equations, we obtain 
rj = - [(1 - b) ln(l - b) + (1 + 6) ln(l + b)} 

(8.11) 



- 2a In 

Solving 7 in (14. 2 p also gives 



(8.12) 



+ a In 



(1 


-6) ln(l 


-6) + (l- 


f 6) ln(l - 


+•6) 






b 2 






(1 


-b) ln(l 


~b) 






(1- 


-b) ln(l - 


-6) + (H 


- 6) ln(l 4 


-b) 



+ 0(a 2 



b 2 



+ 0(a 2 ). 



A combination of (EH}, f lSlTD . (18121) and (E7J) yields 



-l) n+1 r(n + N + 2)N x n- 2x ~ 2 T(x + 1) 



(8.13) 

cf. jnu). 



x < sin 7tx 



r(A^ + l)vr 
1 



1 + 



AT 



+ 0(e 
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9 REGIONS OF VALIDITY FOR EXPANSIONS ( ED ) 
AND ( El ) 



We first consider the case of expansion (18.11) . which is valid for a < 0. Similarity of the two 
asymptotic formulas in (18. 6 p and (I8.13P suggests that the region of validity of expansion 
(18. ip may be extended beyond the origin. This is indeed the case, and we shall show that 
expansion (18. ip actually holds for a G (— oo, a_ — 5] for any small number 5 > 0, a_ being 
the point on the positive real line given in H3 .31) . where the two saddle points w± coalesce 
with each other. 



To prove this, we recall that when < a < a_ — 6, we have f(t (w + ), w + ) < 
f(to(w-), see a statement following (I3.10p . In other words, in the case < a < a_ —5, 
the contribution from the saddle point w + is exponentially small in comparison with that 
from the saddle point So, up to an exponentially small term (E.S.T), we need deal 
with only one saddle point, namely, By an argument similar to that in case (iii) of £J7] 
(see, especially, equation (17.81) ). the integral representation for t n (x, N + 1) given in (13.81) 
can now be written as 



Ux,N + l) 



(9.1) 

where 
(9.2) 



;-i) n r(n + iv + 2) 

7T T(n + 1)T(N - n + 1) 



+oo 



sin(aiV7r) 



_ e Nf(t ( W l W ) dw + E _g_ T 



W 



dt 
d7 



e-^ 2 dr, 



f(t (w), w) =Mn(l - t (w)) + (l-b) IntoH 

+ alnu> — aln(l — w) + b\n[l — (1 — to(w))w], 



which is real on the integration path of w. [Note that the term aln(w — 1) in f(to(w),w) 
is replaced by aln(l — w) in f(to(w), w).] Similar to (I3.27p . we make the transformation 
w — > u defined by 



(9.3) 
with 
(9.4) 

Inserting (19 .3p into (19. ip gives 



f(t Q (w), w) = a ln(u) -u + j, 



a. 



t n {x,N + l) 



(9.5) 

where 
(9.6) 



-l) n T(n + N + 2) 

n r(n + i)r(jv-n + i; 



sin(aA^7r)e 



iV-) 



h(u, r) 



u 



h(u, t) 



e N(aHu)-u) du + E _g_ T 



u dw dt 



-^ 2 dr, 



w — 1 du dr 
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The derivative dw/du is again given by (I3.3ip for u ^ a and (I3.32p for u = a. Let 
h (u,r) = h(u,r), and define the sequence {hi(u, r)} as in fl6.20bj) . Following the proce- 
dure as in case (ii) of 36l we arrive at 



t n (x,N + l) ~ 

(9.7) 



(-l) n r(n + A^ + 2) 
T(n + l)T(N-n + 1) 

' N- aN e N ^T(aN) sin(aA%) y> ' 

Z=0 iV 



where the coefficients q are as given in (16.221) . Although the constant 7 in (19. 7p is 
determined by the mapping defined by (19.31) and (19.41) while the corresponding constant 
7 in (18. ip is determined by a different mapping defined by (I3.27P and (13.281) . the two 
constants are actually the same and can be defined by a single equation 

(9.8) 7 = Re (f(t (w-),w-) - aln(-a)) - a, 

where f(to(w),w) is the same function given in (13.271) or (12.101) . [Recall that 7 is a real 
number by Theorem [1] in §U] Thus, the two expansions (18.11) and (19.71) are exactly the 
same, except for an exponentially small term. The region of validity for expansion (18.11) 
is therefore extended from (—00, 0) to (—00, a_ — 5], < 8 < a_. 

Next, we consider the case of expansion (18. 2p . In the derivation of expansion (18.1 ft when 
a < 0, we have used a mapping w — >• u defined by (I3.27p . The right-hand side of this 
equation is the phase function of the Hankel loop integral for Gamma function; see (I3.25p . 
But, as an alternative, we can also use the phase function in the integral representation 
of the confluent hypergeometric function, as in what we have done when < a < |. This 
will enable us to extend the validity of expansion (18. 2 p to negative values of a. To prove 
this, we mention that for Rec >Re<i, instead of (13. lip , we have 

(9.9) M(d, c, z) = - r (c) ^ (1 ~^ / (-u) d ~\l - uy-^e^du, 

2ixiT{c-d) J 12 

where 72 is the integration path given in (12.81) . This integral representation can be estab- 
lished in the same manner as (12. 8p or (12. 5p . Thus, with d = aN + 1 (a < 0), c = 1 and 
z = r/N, M(<i, c, z) := M(d, c, z)/T(c) can be written as 

1 f N${u) 



(9.10) M(aN+ 1,1, vN) = / e 

2m L 2 



72 u ~ 1 ' 



where ip(u) = aln(— u) — aln(l — u) + rju. The saddle points u± of ip(u) are the same as 
those of ip(u) given in (I3.14p 

, . r]± y/ri 2 + 4ari 

(9.11) u ± = . 

It can be shown that the movements of u± when a < and rj < are like those of the 
saddle points w T of f(t (w),w) when a < 0. Also, we have 

(9.12) M'(aN + 1, 1, V N) = — [ ^— dw, 

2ixi J ~ u 1 
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where the derivative is taken with respect to the third variable in M(d, c, z). 

Now, define the mapping w — > u by 

(9.13) f(t (w), w) = ip(u) + 7 = a ln(— u) — a ln(l — u) + r/u + 7 
with 

(9.14) u(w { + ] ) = u { ^\ u{w-) = u+, 

where r\ and 7 are real numbers. (Note that, similar to (I3.17al) and f!3.17bj) . since w + is 



real and w + > 1, we have w + on both upper and lower edges of the branch cut [0, +00) 
along the real line. We denote these points by 10+ , respectively. The situation is the same 
with u^). Coupling (I9.13P and (I3.23p . the double integral representation for t n (x, N + 1) 
takes the canonical form 

t (xN + 1) J~ 1)n nn + N + 2)e^ 
n[ ' + ' 2m T(n+l)T(N-n + i) 

x / / ^}hll e Nip(u)-N^ dud 



- x •/ - „, " 1 



instead of the one given in (13.291) . In (I9.15p . 
(9.16) h(u,r) u - ldwdt 



w — 1 du dr ' 



and 7« is the steepest descent path of ip(u) passing through u+. Let ho(u,r) = h(u,r), 
and define the sequence {h n (u, r)} as in (I6.9p - (l6.1ip . Following the same argument as in 
Case 1 (0 < a < |) of we obtain 



. Ar ^ (-l) n r(n + iV + 2)e^ 

t n (x, N + 1 - K ,\ r(Ar TTY 

1 (n + 1)1 (N — n + 1) 

(9.17) 



M(aN+l,l }V N)J2^rT 

T~n iV 2 



/=0 

OO 



di 



+M'(aN+l,l,riN)J2T±i 



where c\ and d\ are as given in (16. 131) and (I6.14p . respectively. 

Although the constants rj and 7 in (I9.17P and (18. 2p are determined by two different 
sets of equations, namely, (19.1 3ft - f l9~T4|) and (13. 15[) - (13. 16j) . respectively, one can use a single 
system of equations to determine these two constants for a 6 (—00, 00). Indeed, we can 
define them by 

Re f(t (w-),w-) = Re ip(u + ) +7, 
Re f(t (w + ),w + ) = Re ip(u-) +7, 



(9.18) 



where f(t (w),w) and i/j(u) are the same as those given in (I3.15p . (If u>_ and u + , or 
w + and u_, are on the branch cut, we should choose the correspoinding saddle points 
either both on the upper edge or both on the lower edge of the real line in the equa- 
tions of (I9.18P ). Thus, (I9.17P and (18. 2p are exactly the same, including the coefficients in 
these expansions, and the region of validity of expansion (18. 2 p has been extended from 
< a < \ to -00 < a < \. 
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10 SMALL AND LARGE ZEROS 



It is well-known that from an asymptotic approximation of a function, one can deduce an 
asymptotic approximation for its zeros. The easiest way to achieve this is to apply the 
following result given in Hethcote [23J: 

Lemma 2. In the interval [a—p, a+p], suppose f(t) = g(t)+e(t), where f(t) is continuous, 
git) is differentiable, g(a) = 0, m = min \g'(t)\ > 0, and 

E = max \e(t)\ < min {\g(a — p)\, \g(a + p)\} . 

Then, there exists a zero c of f(t) in the interval such that 

E 

c — a < — . 

m 

To state our result, we first recall the numbers a± given in (13. 3p 

1 ± VT^¥ 
o± = , 

where b is a fixed number in (0, 1). Clearly, < a_ < | < a+ < 1. The number b can be 
as close to or 1 as it may be, but it is fixed. Hence, a_ can be close to but does not 
tend to 0. Similarly, a + can be close to 1 but does not approach 1. 



Theorem 4. LetxN+i iS denote the sth zero of the discrete Chebyshev polynomial t n (x, N+ 
1), arranged in ascending order < xn+i,i < xn+i,2 < • ■ • < XAr + i >n _i < xn+i,u- 
If s — 1 < a^N, then we have 

(10.1) x N+ljS = (s - 1) + O j e W(«-)-^(«+)] j 
as N oo, where 

(10.2) ^(w_) - $(u+) = 2aln 1 + ^} + 4 ^ + nJl + Aa/ri. 

1 - y/1 + Aafq 

If s > n — (1 — a + )N , then we have 

%N+l,s = N — XN+l,n-s+l 

= N — n + s + O j>WK«-)-^+)] I 

Proof. The proof is quite straightforward. For a < 0, we know from (16.211) that the discrete 
Chebyshev polynomials do not have zeros. For < a < a_, it follows from ( I7.13p . (I7.14p 
that both M and M' have zeros at approximately x = m with a exponentially small 
error, where m is a nonnegative integer; see the remark following (I7.13p . The results now 
are obtained by a direct application of Lemma [2] and the symmetry relation (II. 8p . Also, 
note that for < a < a_, the exponent N\ip{uJ) — ip{u + )} in (110.11) is negative, unless 
yT + 4a/ 7] — > in which case a — > a_ (i.e., x is no longer small). □ 



(10.3) 

as N — >■ oo. 
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